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For a given timelike displacement vector the covariant Hamiltonian quasi-local energy expression 
requires a proper choice of reference spacetime. We propose a program for determining the refer- 
ence by embedding a neighborhood of the two-sphere boundary in the dynamic spacetime into a 
Minkowski reference, so that the two sphere is embedded isometrically, and then extremizing the 
energy to determine the embedding variables. Applying this idea to Schwarzschild spacetime, we 
found that for each given future timelike displacement vector our program gives a unique energy 
value. The static observer measures the maximal energy. Applied to the Friedmann-Lemaitre- 
Robertson- Walker spacetime, we find that the maximum energy value is nonnegative; the associated 
displacement vector is the unit dual mean curvature vector, and the expansion of the two-sphere 
boundary matches that of its reference image. For these spherically symmetric cases the reference 
determined by our program is equivalent to isometrically matching the geometry at the two-sphere 
boundary and taking the displacement vector to be orthogonal to the spacelike constant coordinate 
time hypersurface, like the timelike Killing vector of the Minkowski reference. 

PACS numbers: 04.20.Cv, 04.20.Fy, 98.80.Jk 



I. INTRODUCTION 

It is well-known that how to define gravitational en- 
ergy is still an outstanding problem. Since Einstein first 
published his general theory of relativity investigators 
have put much effort into this problem, and some signifi- 
cant progress has been made. Many proposed expressions 
for a gravitational energy-momentum density [l|-|7|, but, 
caught by the equivalence principle, see, e.g., §20.2 in Q, 
they are all pseudotensors rather than tensors. The mod- 
ern concept, introduced by Penrose @, is that properly 
energy-momentum is quasi-local, being associated with 
a closed surface bounding a region (for a comprehensive 
review see 10]). Various ideas about how to define a 
quasi- local energy have been proposed [H, [TTrl23j . Not 
surprisingly the different definitions generally give dif- 
ferent results [24] . In this work we consider one based 
on the covariant Hamiltonian formalism [25l - |29j . In that 
formalism one needs to choose a displacement vector and 
a reference spacetime to determine the ener gy; although 
there was a little exploration of the issue in [26| , no spe- 
cific algorithm for choosing them was prescribed. 

It is generally accepted that the gravitational en- 
ergy for an asymptotically fiat system should be non- 
negative [13, HH . Motivated by this property, it is rea- 
sonable to presume that a suitable reference should be 
the one which minimizes the energy. The results of our 
earlier explorations of this proposal (along with an alter- 
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native analytic approach) are reported in (32^ . A brief 
letter summarizing our later investigations has already 
appeared [331 ] . The present work is intended to be the 
first in a series giving a detailed comprehensive expo- 
sition of our ideas. As in [33], in this paper we will 
embed a two-sphere boundary S and its neighborhood 
in the dynamic spacetime at a constant time and ra- 
dius into the Minkowski reference — with the 2-geometry 
embedded isometrically — and then extremize the energy 
to fix the embedding, and thereby determine the refer- 
ence. We use the Schwarzschild spacetime to test our 
program. A number of investigators have used differ- 
ent quasi-local energy definitions to calculate the energy 
of Schwarzschild spacetime, many of them obtained the 
same results outside the horizon. However, there is no 
such consensus inside the horizon [33TI371 ] . We should 
mention in particular that in [36| the authors, using the 
Brown- York expression, not only get exactly the same re- 
sults for radial geodesic observers as we do but also have 
the same physical interpretations about the observer de- 
pendence and the negative energy result. We show here 
that when applying this energy-extremization program 
to Schwarzschild spacetime for each future timelike dis- 
placement vector there is a uniquely corresponding en- 
ergy but an equivalence class of the references. The ref- 
erence is uniquely determined if we also impose the nor- 
malization condition of the displacement vector in the 
reference spacetime. 

Our program allows us to calculate the energy mea- 
sured by any physical observer. For measuring energy 
inside the black hole horizon, we first show the results 
measured by a radial geodesic observer. Then we exam- 
ine several physical observers, we find they all get rea- 
sonable but different results inside the horizon. If we 
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further vary the energy with respect to the displacement 
vector, then we get the maximum energy, which turns out 
to be that measured by the static observer. To further 
support our program we then consider a slightly differ- 
ent approach. We assume that the displacement vector 
is the timelike Killing vector in the Minkowski reference, 
and then do the extremization. This slightly different 
program produces the same maximum energy as that 
produced by the original program. Although the results 
seem very reasonable, one might still wonder what does 
this energy-extremization program really mean geomet- 
rically. We show that it is equivalent to matching the 
geometry near the two-sphere boundary, i.e., the refer- 
ence is 4D isometric with the dynamic geometry on the 
boundary and it makes the displacement vector N the 
unit normal of the constant coordinate time hypersur- 
face, just like the timelike Killing vector of the Minkowski 
reference. This fact gives us more confidence in the pro- 
gram. 

Applying this program to dynamic spherically sym- 
metric spacetimes, one of the most interesting cases is 
surely homogeneous and isotropic cosmology, i.e., the 
Friedmann-Lemaitre-Robertson- Walker (FLRW) space- 
times. Many researchers have used various pseudoten- 
sor expressions both in general relativity and teleparallel 
gravity to calculate the energy of the FLRW spacetimes 
(see, e.g., [38[ and the references therein). Afshar [3!| cal- 
culated the quasi-local energy of FLRW spacetime by us- 
ing Brown- York's expression and got interesting results. 
In f4(J] the energy value determined for homogeneous cos- 
mologies by the favored covariant Hamiltonian bound- 
ary term with a homogeneous reference was found to be 
zero for all Bianchi class A models and negative for all 
class B models. In this paper we are going to show that 
the energy-extremization program likewise works well in 
dynamic spherically symmetric spacetimes and we give 
some discussion of our results. 

In the following section we briefly review the covari- 
ant Hamiltonian formulation. Then we apply our pro- 
gram to Schwarzschild spacetime in general coordinate 
systems in section 3 and to FLRW cosmology in section 
4; the meaning of the energy extremization program re- 
sults is discussed at the end of each of these sections. The 
conclusion includes a summary of the work. 



II. HAMILTONIAN FORMULATION 

In this section we briefly review the covariant Hamilto- 
nian formalism as developed by our research group 
[29j (for some additional developments along similar lines 
see [4l|, and then apply it to Einstein's gravity the- 
ory: general relativity (GR). 

A first order Lagrangian 4-form for a fc-form field ip 
has the form 



Its variation 

dC dC 

dC = dtdip A p) + dip A — + — A dp, (2) 
dip dp 

defines a pair of first order Euler-Lagrange expressions, 
which are explicitly given by 

^ := dip - d p A, ^ := -<;dp - cLA, (3) 
dp dp 

where <; := (— l) k . The integral of C associates an action 
with any spacetime region. The variation of this action 
is given by the integral of dC. The total differential term 
in ([2]) then leads to an integral over the boundary of 
the region. Hamilton's principle — that the action should 
be extreme with dip vanishing on the boundary — yields 
the field equations: the vanishing of the Euler-Lagrange 
expressions ([3]). 

The action should not depend on the particular way 
points are labeled. Thus it should be invariant under dif- 
feomorphisms, including infinitesimal diffeomorphisms — 
which correspond to a displacement along some vector 
field N. From a gauge theory perspective such displace- 
ments are regarded as "local translations" . Under a local 
translation quantities change according to the Lie deriva- 
tive. Hence, for a diffcomorphism invariant action the 
relation ^ should be identically satisfied when the vari- 
ation operator d is replaced by the Lie derivative £n 
(= diN + *Arrf on the components of form fields): 

dC dC 

d%NC = £^C = d{£Npf\p) + £Npf\~r. — \--=-A£np. (4) 

dip dp 

This simply means that £ is a 4-form which depends on 
position only through the fields p,p. (For this to be the 
case the set of fields in C necessarily includes dynamic 
geometric variables, which means gravity.) 

From the identity ((4]) it directly follows that the 3-form 



H(N) := £ N p Ap-i N C, 



(5) 



satisfies the identity 



8£L 8/L 

- dU(N) = £ N pA — + — A £ N p, (6) 
dtp dp 

showing that it is a conserved "current" on shell (mean- 
ing: when the field equations are satisfied) . Substituting 
the Lagrangian 4-form (JXJ) into the definition ([5]) leads 
to the explicit expression H(N) = d(iN<p A p) + ? i^ip A 
dp + sdtp A ijsip + ijyA, from which one can see that this 
conserved Noether translation current can be written as 
a 3-form linear in the displacement vector plus a total 
differential: 



H(N) =: N^y.f,+dB(N). 



(7) 



Compare the differential of this expression, dH(N) = 
dN^ + N^dHn, with the identity ©; equating the 
dN^ coefficient on both sides reveals that 



C = dp Ap — K{ip,p). 



(1) 



SC dC 

U^-z ll ipA-+,-M ll p 



(8) 
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This (Noether's second theorem type) identity is a neces- 
sary consequence of local diffeomorphism invariance (i.e., 
a symmetry for non-constant 2V*). From this relation one 
can see that vanishes on shell; hence the value of the 
conserved quantity associated with a local displacement 
N and a spatial region £ is determined by a 2-surface 
integral over the boundary of the region: 

E(N,T,):= [ H{N)=(f B(N). (9) 

The value is quasi-local (depending only on the values of 
the fields on the boundary). For any choice of N this 
expression defines a conserved quasi-local quantity. 

One can adjust by hand the 2-form B(N) without af- 
fecting the above argument. Of course this will mod- 
ify the value of the quasi-local quantities. However this 
freedom in choosing the boundary term is entirely under 
physical control, since H(N) is not merely the transla- 
tional Noether current 3-form but is also the Hamiltonian 
3- form which generates the changes (given by £ n) m the 
dynamical quantities. The boundary term in the varia- 
tion of the Hamiltonian reflects the necessary boundary 
conditions implicit in the Hamiltonian. Consequently 
changing B(N) will likewise modify the boundary condi- 
tions implicit in H(N). Thus different choices of bound- 
ary term are associated with different physical boundary 
conditions. 

The above applies quite generally to any geometric 
(i.e., diffcomorphically invariant) gravity theory. Here 
we confine our attention specifically to Einstein's gravity 
theory, general relativity (GR) . GR can be formulated in 
several ways. For our purposes the most convenient is 
to take the orthonormal coframe t9 M = -d^kdx k and the 
connection one-form coefficients oj a p = T a pkdx k as our 
geometric potentials. Moreover we take the connection to 
be a priori metric compatible: Dg a p := dg a p —u>^ a g 1 p — 

pg a ~f = 0. Restricted to orthonormal frames where the 
metric coefficients are constant, this condition reduces to 
the algebraic constraint u; Q/3 = ui^K 

We consider the vacuum (source free) case for simplic- 
ity (the inclusion of sources is straightforward) . GR can 
be obtained from the first order Lagrangian 4-form 

£ GR = ir xP Ap a p+D^AT„-V aP A(j> afi ~t kl p), (10) 

where £l a p := dw a p +u; Q 7 Aw 7 ^ is the curvature 2-form, 
W := d-d^ + uj^ u A 3" is the torsion 2-form, and we 
have made use of the convenient dual form basis rf 1 ^--- ■= 
*{<d a A A • • •)• The 2-forms ft Q/3 , V af} and p af} are 
antisymmetric. We take k :— 8ttG/c 4 = 8ir. Several 
possible boundary terms were identified, each associated 
with a distinct type of boundary condition. In 29] a 
"preferred boundary term" (it has a certain covariant 
property, directly gives the Bondi energy flux, and has a 
positive total energy proof) for GR was identified: 

B(N) = -^(Ar Q ,3 A l nV J + D p N a Ar,J), (11) 



where A indicates the difference between the dynamic 
and reference values, and Dp is the covariant derivative 
using the reference connection. The reference values can 
be determined by pullback from an embedding of a neigh- 
borhood of the boundary into a suitable reference space. 
Now we can use this expression to calculate the gravita- 
tional energy in general relativity. 



III. STATIC SPACETIME: SCHWARZSCHILD 
SPACETIME 

A. Energy-Extremization Program 

Our objective is to find a principle for determining the 
displacement vector and reference spacetime so that we 
can calculate the quasi-local energy for gravitating sys- 
tems. To achieve that we have proposed a program in- 
cluding isometric embedding and energy extremization. 
Now we're going to use the Schwarzschild metric to test 
this program. 

The Schwarzschild metric in the standard spherical co- 
ordinate form is given by 

ds 2 = -Adt 2 + A~ x dr 2 + r 2 dnj, (12) 

where A = l-2ro/r and dfl 2 , = d9 2 +sin 2 Qdtj) 2 - However, 
there are several other reasonable coordinate choices for 
the Schwarzschild metric, such as Painleve-Gullstrand, 
Eddington-Finkelstcin, and Kruskal-Szekeres. In order to 
accommodate most well-known coordinates, we consider 
a more general version of the Schwarzschild metric via a 
coordinate transformation t — t(u, v),r = r(u, v) 

ds 2 = ~(At 2 u - A~ x r 2 u )du 2 + 2(A~ 1 r u r v - At u t v )dudv 
+ (A- 1 r 2 v - At 2 v )dv 2 +r 2 dfl 2 2 . (13) 

Note that under this coordinate transformation we have 
assumed that the orientation is preserved so that the Ja- 
cobian is positive, i.e., \/—a := t u r v — t v r u > 0. Choose 
Minkowski spacetime as the reference: 

ds 2 = -dT 2 + dR 2 + R 2 d0 2 + R 2 sin 2 6d$ 2 . (14) 

The essential issue of the reference choice is the identifi- 
cation between the reference and physical spacetime co- 
ordinates. A legitimate approach for the spherically sym- 
metric case is to assume T = T(u, v), R = R(u, v), O = 
9, $ = (p along with i?o := R(to, ro) — ro; this symmetri- 
cally embeds a neighborhood of the two-sphere boundary 
S at (to , ro ) into the Minkowski reference such that the 
two-sphere boundary is embedded isometrically. 
For any given future timelike displacement vector 

N = N U 3 U + N V 8 V = N l d t + N r 8 r = N T 8 T + N R d n , 

(15) 

We expect that the second term of the expression pT|) 
wont't contribute because the spacetime is spherically 
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symmetric. In particular, if we assume 
3t 



= D T N R = d T N H + T R TT N T + T R RT N R 



= DrN 1 



d R N T 



t T TR N T 



f T RR N R 



d T N 



d B N T , 



(16) 



where all the corresponding connection terms vanish, 
then so does the second term of (fTTj) . Now from (0 fTTj) 
we get for the quasi-local energy associated with a sphere 
of radius r 



E(r) = 



1 

16^ 



(Ar Q ,3 A i NV J + D p N a Ar]J) 



- (N U B + N V C) ■ 



(17) 



where 



B = XT U + g vu (R u ~ 2r u ) + g vv (R v - 2r v ), 
C = XT V + g uu (2r u ~ R u ) + g uv (2r v - R v ), 



(18) 



and the subscripts indicate the related partial differenti- 
ations. Note that the quasi-local energy is evaluated on 
the boundary, the two-sphere S; all the variables appear- 
ing in (fT7"l) and in the following calculations are also eval- 
uated on the two-sphere. Each choice of {T u , T V ,R U ,R V }, 
which we call the embedding variables, means a different 
embedding, hence a different reference. Now extremizing 
the energy with respect to the embedding variables we 
get the conditions 



dE 
Oil 
dE 
df v 
dE 
~dR~ u 

dE 
d~R~ v 



-X 2 T v {N u R u + N v R v ) = 1 

X 2 T u (N u R u + N v R v ) = 0, 

X 2 T V (N U T U + N V T V ) 
+N u g llv - N v g uu = 0, 
-X 2 T U {N U T U + N V T V ) 
+N u g vv - N v g uv = 0. 



(19) 
(20) 

(21) 
(22) 



Note that conditions (|T9l [20)) are equivalent (since we do 
not want both T u and T v to vanish), so we only have 
three independent equations. Using the relations 



g uu = a l g vv , g uv = ~a 1 g uv , g vv = a 1 g uu , 



(23) 



with a = —{t u r v — t v r u ) 2 , then the conditions 
become 

N U R U + N V R V = N R = 0, (24) 
X 2 T v (N u T u +N v T v )-a- 1 (g uv N u +g vv N v ) = 0, (25) 
X 2 T u (N u T u +N v T v )-a~ 1 (g uu N u +g uv N v )=0. (26) 

From X R,, - (|25|) x R u we get 

X(N U T U + N V T V ) + a^ 1 [(g uv N u + g vv N v )R u 

-{g U uN u + g uv N v )R v ] =0. (27) 



Using condition (|24| we get 



N v 



(28) 



and 



iV T = N U T U + N V T V 
N u 



N u 
R v 



(T U R V T v R u ^j 



X R v 

then R v can be found from J27] 



(29) 



Rv N uyy ' v g(N, N) v ' 

Note that we require the displacement vector to be future 
timelike, i.e., N T > and N u > 0, and the orientation 
to be preserved, i.e., the Jacobian is positive. From con- 
dition (|2T)1) we know that R v should be positive, so 



Rv — 



a 



g(N,N) 



N" 



g(N,N) 



N v . (31) 



Now we can calculate the energy. Using the condi- 
tions ([291 EU) we get 



N U B + N V C 
N u 



R r 



y/ag(N,N) + 2(g uv N u r u 



-9uuN u r v + g vv N v r u - g uv N v r v ) 
= 2(t u r v - i„r u ) _1 (y/~g(N, N) - AN 1 ^ , (32) 

where we have used the metric in the calculation. Choose 
N to be unit timelike on the two-sphere: 



l=g(N,N)\i 



= 9uu{N u Y + 2g uv N u N v +g vv (N v Y, (33) 

then we get the quasi-local energy for any given future 
timelike displacement vector N: 

E CX (N) = r [1 - A{N u t u + N v t v )} = r (l - AN*) , 

a result which is independent of any u, v coordinate sys- 
tem choice. It is also independent of the two embedding 
variables T U ,T V . Indeed, we cannot solve for all four vari- 
ables, since we only have three independent conditions. 
So this program produces a unique energy but an equiv- 
alence class of references for any physical observer. 

However, it is reasonable to impose the normaliza- 
tion condition of the displacement vector in the reference 
spacetime, i.e., 

-l=g(N,N) = -(N T ) 2 + {N R ) 2 . (35) 

which together with implies N T = 1. Then 

from dHJ) and (31]) we get 



X 



N u 
R,, 



1 



(36) 
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Using this condition along with 
and T v we get 



to solve for T„ 



T u = At u N l - A-^uW, 
T v = At v N l 



A-\ v N r . 



(37) 



In this way we can calculate for any given observer our 
Hamiltonian boundary term quasi-local energy of the 
Schwarzschild spacetime. 

Now if we further vary the energy with respect to the 
displacement vectors, it should tell us which observer 
would measure the extreme energy. So from 



g(N, N) = ~A{N 1 ) 2 + A-^N 1 ') 2 = -1, 
we take 

a/A/V 4 = coshz, — !=iV = sinhz. 

Va 



(38) 



(39) 



Then we find 
dE 

dz 



cosh z — 1 
1 

E max = r ( 1 



N r = 0, 
A' 



and 



d 2 E 



dz 2 



< 0. 



(40) 



(41) 



"FT " 



So among all physical observers the static observer, i.e., 
energy value 



N = -^jdt, measures the maximum energy (|4T))) . This 



E — 



r (l - v 7 ! - 2m/rj = - 



2m 



+ y/1 - 2m/r 



(42) 



is a standard result [2J, [25|, [2J, [3J, [3J, |43j, |4J . 

Now let us instead assume that the displacement vec- 
tor is the timelike Killing vector in the Minkowski refer- 
ence and then extremize the energy. From the previous 
sections we might already expect that the result should 
be also the value (|40p . Starting from our energy expres- 
sion p7|) we assume 



N = d T = N u d u + N v d v , 

N u = XRv, N v = -XRu, 

-l = g(N,N) 

= 9uu(N u f + 2g uv N u N v + g vv (N v ) 2 
= X 2 {g uu R 2 — 2g uv R u R v + g vv R^) 



(43) 



X 1 



-1 



guuR 2 ^guvRuRv gvvRu 

Then the general energy expression (fl~7f becomes 
l-a^X- 2 + 2a~ 1 {g uv r u R 



(44) 



E= -X 
2 



'-a, 



(45) 



where (|43l l44l have been used. Now there are only two 
variables in our energy expression. We extremize the 
energy with respect to these two variables: 



dE 



dE 
~d~R v 



= =*> (g vv R u - g uv Rv)(l + a^X^ 2 ) 

+2R v (r u R v - r v R u ) = 0, 
= (guuRv - guvRu){l + a^X- 2 ) 



(46) 



-2R u [r u Rv-r v R u ) = Q. (47) 

The combination (|4"oT) x R u + (|4T|) x R v implies 

= (l + a^X-^iguuRl^guvRuRy+gwRl) 
= -X~ 2 (l + a^X^ 2 ) (48) 

X = V-a- 1 = (t u r v - t v r u )-\ (49) 
Now from (j46l |47|) we learn that 

r u R v - r v R u = 0. (50) 
Together with (|44|) this gives 



a — —X 2 — g U uR 2 v — 2g uv R u R v + g V vR 2 u 



R v 



^) In 



so we get 



~A' R * = A- 



(51) 



Pick the plus sign, then we get the unique energy pro- 
duced by this program: 



'■( J - v A) ■ .E mi , s : E, x ( ^= ] (52) 



as expected. 



B. Quasi-local energy measured by other physical 
observers 



It is obvious that the energy value (|4T)|) is only valid 
outside the black hole horizon, since there is no static 
observer inside. Since our energy formula Q34p can be ap- 
plied to any physical observer, to measure the energy in- 
side the black hole let us first examine the radial geodesic 
observers in the Schwarzschild spacetime. For an ob- 
server who falls initially with velocity vo from a constant 
r = a, there would be two different orbits available: (1) 
the crash orbit, where ingoing observers crash directly 
into the singularity at r — and outgoing observers who 
shoot out to the turning point r max and then return back 
to crash; (2) the crash/escape orbit, where ingoing ob- 
servers crash and outgoing observers escape to infinity. 



We choose the unit tangent of the geodesic to be the 
suitable displacement vector, i.e., 



N = 



1 - 



1 2m Ai 2 ' 

1 - — v 1 - v o 



(53) 



where 2m < a and < r < r max . The Hamiltonian 
boundary term energy value for this observer is 



E = r\l- 



^ 2m 



(54) 



This result also agrees with that found using the Brown- 
York quasi- local energy expression [36], [37[ • One can see 
that the energy decreases as the initial velocity vq in- 
creases. When the initial velocity vq is less, equal or 
greater than W2rri/a, which is the escape velocity from 
the Newtonian point of view, the energy is positive, zero 
or negative, respectively. 

The negative value for the energy may appear odd, 
but it can be explained physically. Just as was found in 
another instance , it is correlated with the geometric 
property that the scalar curvature of the spacelike hyper- 
surface orthogonal to the displacement vector is (unlike 
the usual cases) negative. 

To see that let us consider a hypersurface orthogonal 
to the displacement vector 



N 



1 



-di 



2m 



2m. 



(55) 



1/2 



Then the induced metric on the hypersurface is 

h lu , = g lu/ + N^N v , (56) 

where g^ is the four dimensional spacetime metric. Us- 
ing this induced metric one can compute the scalar cur- 
vature R of this hypersurface, which turns out to be 



R 



O 2m 2 

r 2 1 — Vq 



(57) 



One can see that the scalar curvature of this hypersurface 
becomes negative when the initial velocity is greater than 
the escape velocity, i.e. 



vo > 



R<0. 



(58) 



Note that the ingoing geodesic observers can measure 
energy inside the black hole. This energy is proportional 
to the radial distance r, so it is a smooth function in the 
region < r < r max . 

In addition to the static and radial geodesic observers 
there are other natural choices — notably the unit normal 
of the constant coordinate time hypersurface in various 
coordinate systems. 



As an example consider the ADM form of (ingoing) 
Eddington-Finkelstein coordinates [32| : 



ds 2 = — 



du 2 



Ira 



2m , 
1 H dv - 

r 



2mdu 



rjl 



(59) 



where 



2m 

du = dt H A~ 1 dr, dv = dr 

r 

2m 1 

—$*tu 1? tv i 0' 1. 

r 

For this choice the time-component of the coframe and 
the displacement vector are 



du 



\ _j_ 2m 



2 m 



N = e = \ 1+ — d u - 



2m 



2m 



9y, 



consequently 



AT* = N u t u + N v t v = A' 1 (\ + ^\ 2 , (60) 

and the energy (|3"4"]l then has the value 
E EF = r (l - (1 + 2m/r)~ 1 / 2 ) 

= 2m(l + 2m/r + y/l + 2m/r) . (61) 



For (ingoing) Painleve-Gullstrand coordinates: 



ds* = -dv? + \ dv + \l — du] +r 2 dQ^, (62) 



where 

/ 2m 

du = dt + A~ 1 J dr, dv = dr. (63) 

The time-component of the coframe and the displacement 
vector are 

2m 



d° = du, N = e Q =d u - \j— d v , 
2m 



=$>t u = l, t v = -A~ 1 J — , r u = 0, r v = l. 
consequently 

AT* = iV% + N v t v = A -1 , (64) 
and the value of the energy (f34|) now works out to be 

Epo = 0. (65) 
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For this special case the energy vanishes for all r (note 
that the ADM total energy also vanishes for this form of 
the metric). 

In Kruskal-Szekeres coordinates 



ds 2 = e-^(-aiu 2 +aiv 2 )+ r 2 dft 2 , 



N = 



e 4 ™<9„, 



32m 3 

where for r > 2m, u, v are 
r 



(66) 



u = | 1 e 4m sinh - — , 

v 2m / 4m 

r - | l) 2 e~ cosh — — , 

^ 2m / 4m 



(67) 



and we find N l = A 1 / 2 cosh T ^, so the value of the 
energy (|3"4")l is 



£'ks = /■ I I - 1 — cosh 

r 4m 



(68) 



whereas for r < 2m, u, w are 



/ r \ 2 _r_ t 

u = 1 e 4m cosh 



2m J 

— V 

2m/ 



w = I 1 | 2 e 4m sinh 



4m 
t 

4m ' 



(69) 



then N l = {—A) 1 ' 2 cosh ^- and the value of the en- 
ergy ((34]) is 



2m 



t 



E K p, = r I 1 + \l 1 cosh - — 

4m 



(70) 



The expressions for the two regions can be combined to 
give N l = lA]' 1 / 2 cosh^ and 



Eks = r 1 - sgn(A)i/|^| cosh 



4m 



(71) 



These KS results can be compared with those obtained 
by others using the Brown- York quasi-local energy ex- 
pression. They agree with those in (3(| (aside from an 
essential missing absolute value sign) and also (at t = 0) 
with those in |35| . 

Whereas our first results (|42I52|) applied only to the 
region r > 2m where there are static observers, it should 
be noted that all of our freely falling, EF, PG, and KS 
outcomes are well defined, smooth functions of r through 
the horizon. 



C. Observer Adapted Coordinates 

In this section for any coordinate system we take the 
unit normal of the constant coordinate time hypersurface 
as the natural choice of the displacement vector. Now 



using a similar analysis we can find the natural coordi- 
nate system for any observer, which we call the observer 
adapted coordinates. For any N = N t dt + N r d r we can 
find its orthogonal hypersurface. By choosing 

■0° = a dt + aidr, d 1 = b dt 

and from i?°(JV) = l,^ 1 ^) 

2m s 
r 



+ hdr, (72) 



«o 



1 



N\ at 



1 



= we get 
N r , 



-N r : bi = N f . 

(73) 

Assume ?9° = F(t, r)du, where F 1 is an integrating fac- 
tor, then 



dt 



Fdu 



(1 - 2hi) N* 



2m\ 



N r 

w 



-dr. 



and 

ds 2 = -F 2 du 2 - 



1 



1 , N r ' 
— -dr -Fdu 



(74) 



2 dni 



Now for any N we can find the associated observer 
adapted coordinates in which the foliation is along the 
coordinate time, and N is orthogonal to the constant 
time hypersurface. As an example let us consider the 
ingoing geodesic observers: 



N* 

N r 
F = 1. 





1 


>A- 


2m 
a 


1 - 


_ 2m 

r 


vA- 






T - 


2m 
a 







1 - 


v 2 o 



2m 
r 



(76) 



Then 



ds = —du 



dr 



1 - 



^ 2 m 



r 2 dn 2 . 



(77) 



For vq = y/2m/a, which is equivalent to freely falling 
from infinity, this becomes the Painleve-Gullstrand rep- 
resentation. 

There are several different expressions which can 
give m as the result of the quasi-local energy for the 
Schwarzschild spacetime [13, |37[ , so the coordinate sys- 
tem that gives this value should be of interest. By as- 
signing a value m, from our energy expression (134)) we 
find 



N* = (1 - m/r)(l - 2m/r)" 



N r = -m/r, (78) 



and then calculating the integrating factor we get F = 
1 —m/r. From (|75|) we find the form of the metric: 

ds 2 = -(l-m/r) 2 du 2 



(1 — m/r) 1 dr-\ du 



r 2 dVt\. (79) 
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This is the same result as that found in [37] when the 
authors were searching for a foliation such that they can 
identify the Brown- York and the Misner-Sharp energy. 



IV. DYNAMIC SPACETIME: FLRW 
COSMOLOGY 

A. Energy-Extremization Program 



D. The Meaning of Energy Extremization 

As an alternative to extremizing the energy, using (|13[) 
we can instead require that the dynamic and reference 
frames completely match on the 2-surface: 

N = N u d u + N v d v = e = d T , 
i?° = aidu = t?° = ip*dT, 

= aidu + a^dv = fl 1 = ip*dR, 



(80) 



here ip* means the pullback of the embedding map. For 
this kind of isometric matching requirement we find the 
conditions 



T u -- 
Ru - 



Ol, 

= aT 



T v = 0, 
R v — a 3 . 

N v = - 



ai03 



(81) 



From -(tf ) 2 + (d 1 ) 2 + r 2 dSl\ = ds 2 we then get 

{A- l r 2 v - At 2 v )- l / 2 {t u r v - t v r u ), 
{A- X r u r v - At u t v ){A- x rl - At*)" 1 / 2 , 



«i 

«3 



and 



{A-\l-Atlf/ 2 . 



a 2 



2 ■ 

a(a 3 



a, a 



(82) 



(83) 



1^3 



Inserting these expressions along with the results (|8~TI 
52]) into the energy expression (fTT)) gives the isometric 
matching energy: 



E isn = r 



l-rM-'rl-Atly" 2 



(84) 



We now note that this energy value, with further use 
of (I8T1 [82l can be rewritten as 





1 - 




a 2 t v \ 1 








aia 3 J 


= r[ 


1 - 


(N u t u H 


- N v t v )A] 


= r(l - 


N*A) = 





(85) 

i.e., it turns out to have exactly the same value as the 
energy (|34p found from the extremization program. 

Now we can understand what the energy-extremization 
program has done for us: for any unit timelike displace- 
ment vector N the program gives us a 4D isometry on 
the two-sphere boundary between the dynamical and ref- 
erence space, and N is the unit normal of the constant 
coordinate time hypersurface, which is very much like the 
timelike Killing vector of the Minkowski reference being 
orthogonal to the constant T hypersurface. 



With the reference determined following the energy- 
extremization program described in the previous section 
and in [33[, we now apply the procedure to the dynamic 
homogeneous isotropic Friedmann-Lemaitre-Robertson- 
Walker (FLRW) spacetime. The FLRW metric is taken 
in the spherically symmetric form, 



ds 2 



-dt 2 + A 2 dr 2 +a 2 (t) r 2 dfl 2 , 



(86) 



where A = a(t)/y/l - kr 2 and dfl 2 = d9 2 + sin 2 9d<j) 2 . 
Consider a more general version of the FLRW metric ob- 
tained by a coordinate transformation t — t(u,v),r — 
r(u, v) then 



ds 2 = -(t 2 u - A 2 r 2 u ) du 2 + 2(A 2 r u r v - t u t v ) dudv 



+ {A 2 r 2 v -t 2 v )dv z + a z r^dn 



2„2. 



(87) 



Note that we assume that the orientation is preserved 
under coordinate transformations so that the Jacobian is 
positive, i.e., t u r v — t v r u > 0. Choose Minkowski space- 
time as the reference 



ds 2 



-dT 2 + dR 2 



R 2 d<d 2 



R 2 sin 2 6d$ 2 . 



(88) 



We assume, in view of the spherical symmetry, T — 
T(u, v),R = R{u, v), 6 = 0, $ = <f). This embeds a neigh- 
borhood of the boundary two-sphere S at (ioj^o) into 
Minkowksi space. With the restriction Ro := R(to,ro) — 
&(to)ro = aoro the embedding is isometric on the bound- 
ary two-sphere. Assume p^|) . so that the second term of 
the quasi-local expression (fTTj) vanishes. For any given 
future timelike vector (fTS)) from the expressions ([HI fTTj) 
we find the quasi-local energy for a sphere of radius r at 
time t to be 



E(r,t) = -(N U B 



N v Cy 



(89) 



where 



D 
C 
X 

a 

f 



XT U + g uv (R u -2f u )+ g vv (R v - 2f ), 
XT V + g uu (2f u - R u ) + g uv (2f v - R v ), 
(T U R V T v Ry/j , 

QuuSvv 9uv ~~ A- i^u^v ^v^u) ? 

ar. (90) 



Extremizing the energy with respect to the embedding 
variables again leads to the conditions (fl~9Tl22"l) . The con- 
ditions (|T9l |2"0")) are equivalent (since we do not want both 
T u and T v to vanish) , so we only have three independent 
restrictions. Using the relations 



fJ 



vv 1 

9 =a g uu , 



(91) 
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conditions (fT9H22|) become 

N U R U + N V R V = N R = 0, (92) 
X 2 T V {N U T U + N V T V ) - &- l {g uv N u + g vv N v ) = 0, (93) 
X 2 T U {N U T U + N V T V ) - &-\g uu N u + g uv N v ) = 0. (94) 

From ([53 )1 X R,, - ([55] ) X R,, we get 

X(N U T U + N V T V ) + a' 1 \(g U vN u + g vv N v )R u 



-(g U uN u + g uv N v )R v = 0. (95) 



Using (jni]) we find 
Ru = 



N v 



N T = N u T 



Rv, 



N V T V 

so the condition (1951 becomes 



]\ru p 

1 l--a-^g(N,N) = => 



R„ 



u\2 



g(N,NY 



(96) 



(97) 



We require N T > 0, N u > 0, and the Jacobians are 
positive. From (|96| we know R v should be positive, so 



Minkowski reference, i.e., N — dr- Using this condition 
and JSa [Ml BSD to solve for T„ and T„ we find 



T — t N* 
T v — t v N 



- A 2 r u N r 
A 2 r v N r 



(102) 



Now we can calculate the quasi-local energy of the 
FLRW spacetime for any physical observer. One obvious 
choice is N — dt, i.e., the co-moving geodesic observer. 
From (|101|) we get 



E ext (d t ) 



1 - v 7 ! - kr 2 



akr 



1 + VI - kr 2 



(103) 



This is the same as the result found from the Brown- York 
quasi- local energy expression [39| and also in earlier work 
by some of the present authors [45J . One can see this en- 
ergy value is positive, zero, or negative when k = 1, 0, — 1, 
respectively. It is noteworthy that the expression (|103[) 
does not give zero energy for the Milne universe, i.e., 
k = —1, and a = t, which is in fact diffcomorphic to 
Minkowski spacetime. 

Furthermore, we can vary the energy with respect to 
the displacement vector. Since the displacement vector 
is unit timelike, g(N, N) = —(N 1 ) 2 + A 2 (N r ) 2 
take 



T, we 



Rv — 



g(N,N) 



N L 



Ru — —\ 



g(N,N) 



N v . (98) 



Now we can calculate the quasi-local energy. Using 
BHD we find 



N U B + N V C 
N v 

Rv 



a" 1 



y/&gjN, N) + 2{g uv N u f u 
N u r v +g vv N v r u -g uv N v r v ) 



= 2 



g(N,N) 



where we have used the metric and 



f u = a u r + ar u = at u r + ar u , 
r v = a v r + ar v = at v r + ar v , 



(100) 



where d = da/dt. Choose N to be unit timelike on the 
two-sphere (|3"3")l . then we get the quasi-local energy for a 
sphere of radius r at time t for any given future timelike 
displacement vector N: 



E cxt (N) = ar ( 1 - v 7 ! - kr 2 N l - 



VI - kr 2 ' 



-.N r 



. (101) 

This result is independent of the coordinate choice and 
the embedding variables T u and T v . Impose the nor- 
malization condition on the displacement vector in the 
reference spacetime g{N, N) — — 1, then we can uniquely 
determine the reference. Since N R = this condition 
means that N is the timelike Killing vector field in the 



JVr = cosh ; 



AN r = sinh 



z, 



(104) 



then we find 

dz 



= => cosh ; 



VT— kr 2 



Vl kr 2 — a 2 r 2 



N 



Vl — kr 2 



VT- 

or 



kr 2 — 

VT 



-A 



kr 2 



a Vl — kr 2 — a?r 
and the extreme energy value is indeed maximum 



d r , (105) 



Eynax OjT 



1 — \/\ — kr 2 — a? 

ar 3 (k + a 2 ) 
1 + Vl — kr 2 — a 2 r 2 ' 



(106) 



d 2 E 



-ar 



=o 



Vl — kr 2 — a?r 2 < 0. 



One can see that this quasi-local energy is nonnegative 



using the Friedman equation k 
the energy density): 



*7T 

3 



pa (where p is 



E n 



f pa 3 r 3 



1 



1 



^-pa 2 r 2 



> 0, 



(107) 



and it vanishes for the Milne universe. The corresponding 
displacement vector (|105|) is a bit complicated, but it is 
actually just the dual mean curvature vector. The dual 
mean curvature vector V\ is defined as 



V± = fcie - fc ei 



(108) 
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where fco and k\ are the trace of the extrinsic curvature 
of the two-sphere boundary S with respect to e and 
ei respectively. We should also mention here the mean 
curvature vector V is defined as 



V = fciei - k e . 



(109) 



These two vectors are independent of the choice of normal 
frames for S. They depend only on S and constitute a 
set of natural normal vectors for S [46j . 

Now let the displacement vector be the timelike Killing 
vector in the Minkowski reference, and then extremize 
the energy. From the preceding we expect that the result 
should be (|105[ I106[) . If this is indeed so our program 
would be even more satisfactory. Assuming N T = dr 
and -1 = g(N, N) once again leads to (|43I44|) : then ([59]) 
becomes 



E= -X 
2 



2CK {^Quv Tu Quu 



(110) 



Now there are only two variables in our energy expres- 
sion. We extremize the energy with respect to these two 
variables: 



dE 



dR u 

dE 
~dR~ v 







= 



{gwRu - 9uvRv)(l + a 1 X 2 ) 
+2R v (r u R v - f v Ru) = 0, (111) 
(guuRv - 9uvRu){l + a^X- 2 ) 
-2R u (f u R v - f v R u ) = 0. (112) 
Moreover, ([TPTj) x R u + ([TT^]) x R v gives 

= (1 + a' 1 X- 2 )(g uu R 2 v - 2g uv R u R v + g vv R 2 u ) 
= -X^il + a^X- 2 ) 

(113) 



x = V-t 



Now pTTl [TT2"|) tell us 

' It R v - f v R 
and together with (|4"4")l we get 

2 f2 

R„ 



o. 



(114) 



9uu^y 2g uv r u v v -\- g^v^u 



1 — kr 2 — a 2 r 2 



R 



1 — kr l — a z r z 

where (|100[) was used. Pick the plus sign, then we get 
the unique energy produced by this program: 



N 



y/T— kr 2 dr y/T— kr 2 

-A 7 =====d r , 



\/l — kr 2 — d 2 r 2 a \/\ — kr 2 — a 



E{d T ) c *t =ar(l- yj\ - kr 2 - 



(116) 



This result agrees with E max (fTTOl [TU71) . 



B. The Meaning of Energy Extremization and Its 
Outcomes 



Instead of extremizing the energy, one can require that 
the dynamic frame associated with the metric (|87p ex- 
actly matches the reference frame on the boundary: 

d T , 



N = N u d u + N v d v = e 
= a x du = = <p*dT, 
•d 1 = a 2 du + d^dv — t? 1 = ip*dR, 



(117) 



where tp* means pullback. From this kind of matching 
we get 



oi, T v 



= 
N v = 



Ru = a>2, 
did 3 



R, 



Ct3, 



(118) 



From — ($ ) 2 + (i} 1 ) 2 + r 2 dfl 2 — ds 2 , and picking the plus 
sign we get 

A(A 2 r 2 v ^t 2 v )-^ 2 (t u r v ~t v r u ), 



d 2 = (A 2 r u r v - t u t v )(A 2 r 2 v - t 2 v )- 1/2 , 

a 3 = (i\ 2 -a i/2 . 

Putting (fTTgl rTT9|) into flM]) then gives 



ar 



l-(A 2 r 2 v -t 2 v )'^ 2 (ar v +drt v ) 



(119) 



(120) 



On the other hand, transforming the value we found 
above in (I101[) yields 



-Eext(#r) = ar 



1 - \/l - kr 2 N l 



s/T— kr 2 



-.N r 



1 - yJl-kr 2 {N u t u + N v t v ) 



(N u r u + N v r v ) 



vl — kr 2 



ar 



l-{A 2 r 2 v -t 2 v )- 1/2 {ar v +drt v ) 



(121) 



So for any unit timelike displacement vector N this pro- 
gram gives us a 4D isometry at the two-sphere boundary 
between the dynamic geometry and the reference, and 
N is the unit normal of the constant coordinate time 
hypersurface — which is similar to the timelike Killing 
vector of the Minkowski reference being orthogonal to 
the constant T hypersurface. 

From (|108l I117[) for this case we find 



fco = 0, k\ = 1. 
And from ./V = 8t = eo we also know 
k = 0, ki = 1, 



(122) 



(123) 



where bar represents the reference quantities. Physically 
this means that the expansions of the two-sphere bound- 
ary S and its image f(S) are the same, since the mean 
curvature vectors (aka expansion vectors) are the same. 
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Many investigators believe that quasi-local energy 
should be positive (e.g., [H, HtJ), but one can see that 
for k — — 1 the quasi-local energy (|103|) is negative. On 
the other hand, the quasi-local value (| 107[) is always pos- 
itive. Thus the quasi-local energy value depends on the 
observer, but there is a set of observers in FLRW space- 
time who would measure the maximum energy — which is 
nonnegative. Furthermore, if we wisely choose the ref- 
erence for an observer, it should be possible to always 
get nonnegative energy values, since the physically mean- 
ingful energy is measured relative to the ground state. 
Although positivity is a nice property for gravitational 
energies, negative energies can also be reasonable, when 
due consideration is given to the dynamics, the differ- 
ent choices of observer and the ground state (for some 
discussion of this topic see Ref. |4G|V 



V. CONCLUSION 

The covariant Hamiltonian quasi-local energy expres- 
sion has several virtues, however it also suffers from two 
ambiguities: which displacement vector and which refer- 
ence. We propose embedding the two-sphere boundary 
and its neighborhood in the dynamic spacetime into the 
Minkowski reference, with the 2-geometry being embed- 
ded isometrically, and then extremizing the energy to fix 
the embedding and thereby the reference. In this paper 
we use the Schwarzschild spacetime and FLRW cosmol- 
ogy to test this idea and obtained encouraging results. 

For each future timelike displacement vector the pro- 
gram produces a uniquely corresponding energy, so that 
we can discuss the energy measured by different ob- 
servers. For example, a radial geodesic observer in the 
Schwarzschild spacetime can measure the energy at any 



distance from the singularity — even inside the black hole. 
The energy measured can be positive, zero or even nega- 
tive. We found that the negative result is related to the 
negative scalar curvature of the 3-hypersurfacc. By im- 
posing some further conditions on the reference we can 
get a positive energy. This program can also be applied 
to the Reissner-Norstrom spacetime or other static spher- 
ically symmetric spacetimes. 

For the dynamic FLRW cosmology, the sign of the en- 
ergy of the co-moving observer is k, the sign of the spa- 
tial curvature, which is negative for the open universe. 
When we vary the energy with respect to the observers 
we find that the maximum is nonnegative. The physically 
meaningful energy is defined relative to a ground state. 
So in this sense, although positivity is a nice property 
for gravitational energies, negative energies can also be 
reasonable, considering the dynamics, the different ob- 
servers and the choice of the ground state. Moreover, 
we find that this program is actually equivalent to iso- 
metrically matching the 4D geometry at the two-sphere 
boundary and making the displacement vector orthogo- 
nal to the spacelike constant coordinate time hypersur- 
face, just like the timelike Killing vector of the Minkowski 
reference. This matching has a clear geometrical mean- 
ing for FLRW cosmology: it matches the expansion of 
the two-sphere boundary S and that of its reference im- 
age ip(S). 
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